Microarray-based gene expression analysis combined with laser capture microdissection is beneficial in investigating the modes of action of ocular toxicity.

The retina consists of several layers, and drugs can affect the retina and choroid separately. Therefore, investigating the target layers of toxicity can provide useful information pertaining to its modes of action. Herein, we compared gene expression profiles obtained via microarray analyses using samples of target layers collected via laser capture microdissection and samples of the whole globe of the eye of rats treated with N-methyl-N-nitrosourea. Pathway analyses suggested changes in the different pathways between the laser capture microdissection samples and the whole globe samples. Consistent with the histological distribution of glial cells, upregulation of several inflammation-related pathways was noted only in the whole globe samples. Individual gene expression analyses revealed several gene expression changes in the laser capture microdissection samples, such as caspase- and glycolysis-related gene expression changes, which is similar to previous reports regarding N-methyl-N-nitrosourea-treated animals; however, caspase- and glycolysis-related gene expressions did not change or changed unexpectedly in the whole globe samples. Analyses of the laser capture microdissection samples revealed new potential candidate genes involved in the modes of action of N-methyl-N-nitrosourea-induced retinal toxicity. Collectively, our results suggest that specific retinal layers, which may be targeted by specific toxins, are beneficial in identifying genes responsible for drug-induced ocular toxicity.

Establishment of an in vitro cytotoxicity assay platform using primary monkey cardiomyocytes.

To establish an in vitro cytotoxicity assay platform using monkey cardiomyocytes, we isolated primary cardiomyocytes from fetal cynomolgus monkeys at different gestation days (from day 39 to 90) using the trypsin and collagenase digestion method, which was identical to the standard procedure for rat cardiomyocytes. Under these conditions, the primary cells obtained from monkeys at gestation day 63 or earlier showed spontaneous beating, with >80% cells being viable from all fetuses. Transcriptome analysis of the monkey cardiomyocytes indicated that the cells have essential components of cardiac functions, such as myosins, α-actin, cardiac troponins, and calcium-related molecules. The susceptibility to doxorubicin-induced cytotoxicity in monkey cardiomyocytes was comparable to that in rat cardiomyocytes, as evaluated based on intracellular ATP levels. Microarray analysis with Ingenuity Pathway Analysis revealed that doxorubicin predominantly increased the expression of several key genes involved in the endoplasmic reticulum stress pathway in monkey cardiomyocytes than in rat cardiomyocytes. In conclusion, we isolated primary monkey cardiomyocytes that showed similar sensitivity to doxorubicin as compared with rat cardiomyocytes. This in vitro monkey cardiomyocyte assay platform would serve as a powerful tool for the investigation of the interspecies differences in drug-induced cardiotoxicity and its underlying mechanism.

Absolute Quantification of Plasma MicroRNA Levels in Cynomolgus Monkeys, Using Quantitative Real-time Reverse Transcription PCR.

RT-qPCR is one of the most common methods to assess individual target miRNAs. MiRNAs levels are generally measured relative to a reference sample. This approach is appropriate for examining physiological changes in target gene expression levels. However, absolute quantification using better statistical analysis is preferable for a comprehensive assessment of gene expression levels. Absolute quantification is still not in common use. This report describes a protocol for measuring the absolute levels of plasma miRNA, using RT-qPCR with or without pre-amplification. A fixed volume (200 µL) of EDTA-plasma was prepared from the blood collected from the femoral vein of conscious cynomolgus monkeys (n = 50). Total RNA was extracted using commercially available system. Plasma miRNAs were quantified by probe-based RT-qPCR assays which contains miRNA-specific forward/reverse PCR primer and probe. Standard curves for absolute quantification were generated using commercially available synthetic RNA oligonucleotides. A synthetic cel-miR-238 was used as an external control for normalization and quality assessment. The miRNAs that showed quantification cycle (Cq) values above 35 were pre-amplified prior to the qPCR step. Among the 8 miRNAs examined, miR-122, miR-133a, and miR-192 were detectable without pre-amplification, whereas miR-1, miR-206, and miR-499a required pre-amplification because of their low expression levels. MiR-208a and miR-208b were not detectable even after pre-amplification. Sample processing efficiency was evaluated by the Cq values of the spiked cel-miR-238. In this assay method, technical variation was estimated to be less than 3-fold and the lower limit of quantification (LLOQ) was 102 copy/µL, for most of the examined miRNAs. This protocol provides a better estimate of the quantity of plasma miRNAs, and allows quality assessment of corresponding data from different studies. Considering the low number of miRNAs in body fluids, pre-amplification is useful to enhance detection of poorly expressed miRNAs.

Comprehensive Analysis of Circulating microRNA Specific to the Liver, Heart, and Skeletal Muscle of Cynomolgus Monkeys.

Circulating microRNAs (miRNAs) could represent sensitive and specific biomarkers for tissue injury. However, their utility as biomarkers in nonclinical toxicological studies using nonhuman primates is limited by a lack of information on their organ specificity and circulating levels under resting condition of the animals. Herein, liver, heart, and skeletal muscle-specific expression patterns of miRNAs were determined in 27 tissues/organs from male and female monkeys (n =2/sex) by next-generation sequencing (NGS) analysis. This analysis revealed organ-specific miRNAs in the liver (miR-122), heart (miR-208a and miR-499a), and skeletal muscle (miR-206). Next, plasma was collected from conscious-naive male and female cynomolgus monkeys (n = 25/sex) to better understand the expressions of organ-specific circulating miRNAs. The absolute values of circulating miRNAs were quantified using a Taqman microRNA assay. MiR-1, miR-133a, and miR-208b showed preferential expression in the heart and skeletal muscles, whereas miR-192 was abundant in the liver, stomach, small intestine, and kidney. These miRNAs had identical sequences to their human counterparts. Six organ-specific miRNAs (miR-1, miR-122, miR-133a, miR-192, miR-206, and miR-499a) could be evaluated quantitatively by quantitative real-time reverse transcription polymerase chain reaction with or without preamplification. No significant sex differences were noted for these circulating miRNAs. For their circulation levels, miR-133a showed more than 900-fold interindividual variation, whereas miR-122 showed only a 20-fold variation. In conclusion, we profiled circulating organ-specific miRNAs for the liver, heart, and skeletal muscle of cynomolgus monkeys.

A monkey model of acetaminophen-induced hepatotoxicity; phenotypic similarity to human.

Species-specific differences in the hepatotoxicity of acetaminophen (APAP) have been shown. To establish a monkey model of APAP-induced hepatotoxicity, which has not been previously reported, APAP at doses up to 2,000 mg/kg was administered orally to fasting male and female cynomolgus monkeys (n = 3-5/group) pretreated intravenously with or without 300 mg/kg of the glutathione biosynthesis inhibitor, L-buthionine-(S,R)-sulfoximine (BSO). In all the animals, APAP at 2,000 mg/kg with BSO but not without BSO induced hepatotoxicity, which was characterized histopathologically by centrilobular necrosis and vacuolation of hepatocytes. Plasma levels of APAP and its reactive metabolite N-acethyl-p-benzoquinone imine (NAPQI) increased 4 to 7 hr after the APAP treatment. The mean Cmax level of APAP at 2,000 mg/kg with BSO was approximately 200 µg/mL, which was comparable to high-risk cutoff value of the Rumack-Matthew nomogram. Interestingly, plasma alanine aminotransferase (ALT) did not change until 7 hr and increased 24 hr or later after the APAP treatment, indicating that this phenotypic outcome was similar to that in humans. In addition, circulating liver-specific miR-122 and miR-192 levels also increased 24 hr or later compared with ALT, suggesting that circulating miR-122 and miR-192 may serve as potential biomarkers to detect hepatotoxicity in cynomolgus monkeys. These results suggest that the hepatotoxicity induced by APAP in the monkey model shown here was translatable to humans in terms of toxicokinetics and its toxic nature, and this model would be useful to investigate mechanisms of drug-induced liver injury and also potential translational biomarkers in humans.

(+)-Usnic acid-induced myocardial toxicity in rats.

(+)-Usnic acid (UA) has been known to be a strong uncoupler, and mitochondrial and endoplasmic reticulum (ER)-related stresses are suggested to be involved in the mechanism of hepatotoxicity. However, it has not been clarified whether UA causes toxicity in other mitochondria-rich organs such as the heart. We elucidated whether UA induces cardiotoxicity and its mechanism. UA was orally administered to rats for 14 days, and laboratory and histopathological examinations were performed in conjunction with toxicogenomic analysis. As a result, there was no alteration in blood chemistry, whereas cytoplasmic rarefaction of myocardium was observed microscopically. This finding corresponded to the swollen mitochondria observed ultrastructurally. Immunohistochemically, expression of prohibitin, indicating mitochondrial imbalance, increased in the sarcoplasmic area. Toxicogenomic analysis highlighted the upregulation of gene groups consisting of oxidative stress, ER stress, and amino acid limitation. Interestingly, the number of upregulated genes was larger in the amino acid limitation-related gene group than that in other groups, implying that amino acid limitation might be one of the sources of oxidative stress, not only mitochondria and ER-originated stresses. In conclusion, the heart was manifested to be one of the target organs of UA. Mitochondrial imbalance with complex stresses may be involved in the toxic mechanism.

Toxicogenomic investigation on rat testicular toxicity elicited by 1,3-dinitrobenzene.

Rats were treated with a single oral dose of 10, 25 and 50mg/kg of 1,3-dinitrobenzene (DNB), and the testis was subjected to a GeneChip microarray analysis. A total of 186 and 304 gene probe sets were up- and down-regulated, respectively, by the DNB treatment, where spermatocyte death and Sertoli cell vacuolation in testis and increased debris of spermatogenic cell in epididymis were noted. The expression profile for four sets of genes were investigated, whose expressions are reported to localize in specific cell types in the seminiferous epithelium, namely Sertoli cells, spermatogonia plus early spermtocytes, pachytene spermatocytes and round spermatids. The data demonstrated that pachytene spermatocyte-specific genes elicited explicit down-regulation in parallel with the progression of spermatocyte death, while other gene sets did not show characteristic expression changes. In addition, Gene Ontology analysis indicated that genes associated with cell adhesion-related genes were significantly enriched in the up-regulated genes following DNB treatment. Cell adhesion-related genes, namely Cdh2, Ctnna1, Vcl, Zyx, Itgb1, Testin, Lamc3, Pvrl2 and Gsn, showed an increase in microarray and the up-regulation of Cdh2 and Testin were confirmed by real time RT-PCR. The gene expression changes of pachytene spermatocyte-specific genes and cell adhesion-related genes were thought to reflect a decrease in the number of spermatocytes and dysfunction of Sertoli-germ cells adhesion junction, and therefore these genes would be potential genomic biomarkers for assessing DNB-type testicular toxicity.

Sex and circadian modulatory effects on rat liver as assessed by transcriptome analyses.

The present study was designed to fully uncover sex and circadian modulatory effects on rat liver. Hepatic transcriptome analyses were performed at 4 hr intervals of a day-night cycle using young adult male and female rats. Sexually dimorphic genes, which were identified by a cross-sex comparison of time series data, included representative sex-predominant genes such as male- or female-predominant cytochrome P450 subfamilies (Cyp2c11, Cyp2c12, Cyp2c13, and Cyp3a2), sulfotransferases, and glutathione S-transferase Yc2. The identified sexually dimorphic genes were over-represented in the metabolism of retinols, xenobiotics, linoleic acids, or androgen and estrogen, or bile acid biosynthesis. Furthermore, transcription factor targets modeling suggested that transcription factors SP1, hepatocyte nuclear factor 4-alpha (HNF4-alpha), and signal transducer and activator of transcription 5b (STAT5b) serve as core nodes in the regulatory networks. On the other hand, Fourier transform analyses extracted universal circadian-regulated genes in both sexes. The circadian-regulated genes included clock or clock-controlled genes such as aryl hydrocarbon receptor nuclear translocator-like (Arntl), period homolog 2 (Per2), and D site albumin promoter binding protein (Dbp). The extracted cyclic genes were over-represented in major tissue activities, e.g. the urea cycle and the metabolism of amino acids, fatty acids, or glucose, indicating that the major liver functions are under circadian control. The transcription factor targets modeling suggested that transcription factors SP1, HNF4-alpha, and c-Myc proto-oncogene protein (c-MYC) serve as major hubs in the circadian-regulatory gene networks. Interestingly, transcription factors SP1 and HNF4-alpha are likely to orchestrate not only sexually dimorphic, but also circadian-regulated genes even though each criterion was rather mutually exclusive. This suggests the cross-talk between those regulations. Sexual dimorphism is likely to interact with circadian rhythmicity via overlapping gene regulatory networks on rat liver.

Hepatic transcriptome and proteome responses against diethyl maleate-induced glutathione depletion in the rat.

Hepatic transcriptome and proteome responses against glutathione depletion were investigated by Affymetrix GeneChip Microarray and 2-dimensional gel electrophoresis (2D-DIGE), followed by MALDI-TOF-MS analysis and utilizing a glutathione-depleted rat model treated with diethyl maleate (DEM). Hepatic glutathione content decreased to 1.29 µmol/g liver (25.5% compared to control) after DEM treatment, and there were no apparent hepatotoxic signs estimated by blood chemistry examinations. A total of 247 and 213 annotated gene probe sets exhibited greater than twofold up- and down-regulation compared with controls, respectively. The up-regulated gene list contained a number of glutathione depletion-responsive genes reported previously, such as Trib3, Srxn1, Myc, Asns, Igfbp1, Txnrd1, or Hmox1, suggesting that these genes are robust mRNA biomarkers for evaluating hepatic glutathione depletion. In the 2D-DIGE analysis, proteins for a total of 361 spots were identified by MALDI-TOF-MS analysis. Of the identified proteins, 5 and 14 proteins showed up- and down-regulation, respectively. Some proteins exhibited differential expression in the protein level but not in the mRNA level, including L-FABP, MAWDBP, aldo-keto reductase family 1 member A1, catalase and ATP synthase subunit beta, suggesting that these proteins would be potential protein biomarkers for evaluating glutathione depletion. Moreover, up-regulation of FABP1 protein along with up-regulation of PPARα-regulated gene transcripts (i.e., Acot2 and Acot4) is indicative of PPARα activation, which may contribute to hepatocellular protection against glutathione depletion-induced oxidative stress. The up-regulation of L-FABP1 was detected by proteome data but not by transcriptome data, demonstrating the advantage of utilizing transcriptomics and proteomics combination to investigate glutathione depletion-induced molecular dynamics.

Circadian modulation of hepatic transcriptome in transgenic rats expressing human growth hormone.

The secretory profile of growth hormone (GH) is sexually dimorphic in rats. In male transgenic (TG) rats expressing human GH (hGH) that we generated, the circulating levels of both hGH and endogenous GH are flattened with no male-type pulsatility. To elucidate the regulatory role of episodic GH profile on the liver, the hepatic transcriptome of male TG rats at the middle of the light and dark phases was characterized by genome-wide analyses as compared with that of male wild-type (WT) rats. Transcripts commonly up- or down-regulated regardless of the lighting conditions in TG rats were mainly enriched in the metabolism of xenobiotics. In TG rats, the gene expression profile was functionally feminized, verifying that the sexually dimorphic profile of GH rather than genetic sexuality is a stronger sex-determining factor on the hepatic transcriptome. The common transcripts which fluctuated during the day in both TG and WT rats were enriched in circadian rhythm signaling, and physiological rhythmicity was considered to be finely interconnected with liver metabolism via sexually dimorphic GH secretion. In contrast, some genes were differentially regulated in TG rats at only one of two time points measured, and others were fluctuated daily in only one genotype. In particular, some genes involved in the GH signaling pathway were included, suggesting the signal transduction is circadian-modulated depending upon the GH profile. Our transcriptome analyses clarified the regulatory role of episodic GH profile on the liver and strengthen the functional link between sexually dimorphic GH secretion, liver metabolism, and its circadian regulation.

Scoring multiple toxicological endpoints using a toxicogenomic database.

As information regarding microarray data sets and toxicogenomic biomarkers grows rapidly, the process of analyzing data and interpreting the results is increasingly complicated. To facilitate data analysis, a simple expression ratio-based scoring method called the TGP1 score was previously proposed [Kiyosawa, N., Shiwaku, K., Hirode, M., Omura, K., Uehara, T., Shimizu, T., Mizukawa, Y., Miyagishima, T., Ono, A., Nagao, T., Urushidani, T., 2006. Utilization of a one-dimensional score for surveying chemical-induced changes in expression levels of multiple biomarker gene sets using a large-scale toxicogenomics database. J. Toxicol. Sci. 31, 433-448]. Although the TGP1 score has demonstrated its efficacy for rapid comprehension of large-scale toxicogenomic data sets, inclusion of low quality gene expression data in the biomarker gene set produced flaws in the calculated score. To overcome this shortcoming, we tested a new scoring method called the differentially expressed gene score (D-score), where Detection Call as well as signal log ratios generated by MAS5 algorithm on Affymetrix GeneChip data were considered for the calculation. Four prototypical toxicants, namely acetaminophen, phenobarbital, clofibrate and acetamidofluorene, were used for detailed analysis. A toxicogenomics database (TG-GATEs) was utilized as a reference data set. The D-score successfully alleviated the effects of low quality data on the score calculation, and captured the overall direction of expression changes as well as the magnitude of expression change level of a set of genes, highlighting the affected toxicological endpoints elicited by chemical treatment. The D-score will be useful for high-throughput toxicity screening using a toxicogenomic database and biomarkers.

Utilization of a toxicogenomic biomarker for evaluation of chemical-induced glutathione deficiency in rat livers across the GeneChip data of different generations.

Previously, we reported 69 probe sets (GSH probe sets) of RG U34A GeneChip that were useful for the evaluation of chemical-induced glutathione depletion in rat livers. The aim of the present study was to investigate whether these probe sets could be applied to the analysis of RAE 230A GeneChip data. Since a straightforward data comparison of RG U34A and RAE 230A GeneChips could not overcome the generation-dependent discrepancy in signal profiles, we tried two methods to improve the data compatibility between the two GeneChips. First, we re-calculated the signal values by excluding the probes with poor-overlapping sequences between the two GeneChips, but the data compatibility did not improve from the view point of Spearman's and Pearson's correlation coefficients. On the other hand, the PCA result demonstrated that an adjustment of the baseline signal level between the RG U34A and RAE 230A GeneChip data on vehicle-treated rats dramatically improved the data compatibility, suggesting that the GSH probe sets identified from RG U34A GeneChip data can be utilized in RAE 230A GeneChip data as well. Such a baseline adjustment of signal data is an easy and practical way to utilize biomarkers across GeneChip data of different generations.

Evaluation of glutathione deficiency in rat livers by microarray analysis.

Hepatic glutathione content was measured and gene expression data were obtained using an Affymetrix RG U34 array after treatment with tap water containing 20mM l-buthionine (S, R)-sulfoximine (BSO) to male F344 rats for four consecutive days. Both Spearman's and Pearson's correlation coefficients were calculated between the glutathione content and the mRNA content level obtained from the microarray analysis individually. Sixty-nine gene probes, which were statistically significant (Spearman's correlation, P < 0.05) and showed a Pearson's correlation coefficients (Pearson's r) less than -0.8 between mRNA content and hepatic glutathione content, were identified as glutathione deficiency-correlated probes. By comparing the hepatic gene expression profiles between BSO- and butylated hydroxyanisole (BHA)-treated rats, 14 probes of genes that showed an increase in the corresponding gene mRNA levels only after the BSO treatment were thought to be good indicators of glutathione deficiency. A principal component analysis successfully illustrated the time-course of hepatic gene expression after the treatment with acetaminophen, phenobarbital and clofibrate, and the expression profiles were thought to reflect the changes in hepatic glutathione levels. The identified gene probes in the present study would be useful as markers for assessing hepatocellular glutathione deficiency, or oxidative stress level, based on microarray data.

Molecular mechanism investigation of phenobarbital-induced serum cholesterol elevation in rat livers by microarray analysis.

Phenobarbital (PB) increases serum total cholesterol levels in rodents and humans. To investigate the underlying molecular mechanisms, we performed a microarray analysis on liver of rats treated repeatedly with 100 mg/kg PB, and examined the serum blood chemistry. The serum concentration of non-esterified fatty acids was decreased from day 1 to day 14 except for day 7, and that of cholesterol was increased from day 4 to day 14. The serum concentration of total ketone bodies was increased on day 7, and that of triglycerides was decreased on day 14. Transcript content of glycolytic genes was decreased by PB treatments, while that of lipoprotein lipase was continuously increased, suggesting a notion that repetitive PB treatments impaired glycolysis and stimulated lipolysis in the liver. The hypothesis was examined by using a previously reported flux-balance model. The increase in mRNA content of malic enzyme after the PB treatment agreed well with the flux-balance model result, suggesting the validity of our hypothesis. The findings also suggested that there was an abundance of acetyl-CoA and shortage of glycolytic products after the repeated PB treatments. Although ketogenesis would normally occur under such cellular conditions, it was only weakly observed after the repeated PB treatments, presumably owing to a decrease in HMG-CoA synthase mRNA content. On the other hand, the mRNA content of several cholesterogenic genes was slightly induced by PB treatments. Thus, serum chemistry and microarray results suggested that repeated PB treatments induced cholesterogenesis in rat livers, which may have contributed to the elevation of the serum total cholesterol concentration.

Effect of serum cholesterol on the mRNA content of amyloid precursor protein in rat livers.

Genes that showed mRNA content profiles, which correlated with serum concentrations of total cholesterol (T.CHO), were screened from the microarray data of phenobarbital (PB)- or clofibrate (CLO)-treated rat livers, and the correlation was evaluated based on Spearman's correlation coefficient. Many genes involved in the cholesterol or bile acid metabolism were highly correlated such as UDP-glucuronosyltransferase-21, apolipoprotein A-I and cMOAT. The mRNA content of the amyloid precursor protein (APP) showed the 5th highest correlation among the 8799 probes in the Affymetrix Rat Genome U34 Array. In the livers of rats fed a high-cholesterol (1%) diet for 33 days, serum T.CHO levels increased by 4.6-fold, and the hepatic APP mRNA content also increased by 1.9-fold compared to the control group. These data suggest that the hepatic APP mRNA content was affected by serum T.CHO, and that hepatic APP was involved in cholesterol metabolism in rat livers.

Phylogenetic tree facilitates the understanding of gene expression data on drug metabolizing enzymes obtained by microarray analysis.

The gene expression data of drug metabolizing enzymes (DMEs) in male F344 rat livers were examined after treatments with phenobarbital (PB), clofibrate (CPIB), 3-methylcholanthrene (3-MC) or butylated hydroxyanisole (BHA) using an Affymetrix GeneChip system. Nucleotide sequence-based phylogenetic trees combined with a heat map, that presents both quantitative and qualitative data, were created. Most DME gene probes were successfully classified into the corresponding gene families, although a few were not due to the presence of non-coding or promoter region sequences in the target gene. There were also some data discrepancies among probes of the same gene family, indicating the inappropriate design of these probes. With this method, microarray probes with confusing nomenclature and quality differences can be identified. In addition, a good correlation between the gene expression data and protein data was confirmed, indicating the usefulness of this method for the comprehensive monitoring of DME activity in rat livers treated with xenobiotics.